Identification and quantitation of multiple variants in RNA virus genomes.

The goal of the study was to identify and characterize RNA virus variants containing mutations spread over genomic distances >5 kb. As proof of concept, high-quality viral RNA of the Dengue 2 component of Takeda's tetravalent dengue vaccine candidate (TDV-2) was used to develop a reverse transcription-polymerase chain reaction protocol to amplify a ∼5.3 kb cDNA segment that contains the three genetic determinants of TDV-2 attenuation. Unique molecular identifiers were incorporated into each viral cDNA molecule for PacBio library preparation to improve the quantitative precision of the observed variants at the attenuation loci. Following assay optimization, PacBio long-read sequencing was validated with multiple clone-derived TDV-2 revertant variants and four complex revertant mixtures containing various compositions of TDV-2 and revertant viruses. PacBio sequencing analysis correctly identified and quantified variant composition in all tested samples, demonstrating that TDV-2 revertants could be identified and characterized and supporting the use of this method in the differentiation and quantification of complex variants of other RNA viruses. Long-read sequencing can identify complex RNA virus variants containing multiple mutations on a single-genome molecule, which is useful for in-depth genetic stability and revertant detection of live-attenuated viral vaccines, as well as research in virus evolution to reveal mechanisms of immune evasion and host cell adaption.

Annotated genome sequence of a fast-growing diploid clone of red alder (Alnus rubra Bong.).

Red alder (Alnus rubra Bong.) is an ecologically significant and important fast-growing commercial tree species native to western coastal and riparian regions of North America, having highly desirable wood, pigment, and medicinal properties. We have sequenced the genome of a rapidly growing clone. The assembly is nearly complete, containing the full complement of expected genes. This supports our objectives of identifying and studying genes and pathways involved in nitrogen-fixing symbiosis and those related to secondary metabolites that underlie red alder's many interesting defense, pigmentation, and wood quality traits. We established that this clone is most likely diploid and identified a set of SNPs that will have utility in future breeding and selection endeavors, as well as in ongoing population studies. We have added a well-characterized genome to others from the order Fagales. In particular, it improves significantly upon the only other published alder genome sequence, that of Alnus glutinosa. Our work initiated a detailed comparative analysis of members of the order Fagales and established some similarities with previous reports in this clade, suggesting a biased retention of certain gene functions in the vestiges of an ancient genome duplication when compared with more recent tandem duplications.

High-Quality Reference Genome Sequence for the Oomycete Vegetable Pathogen Phytophthora capsici Strain LT1534.

The oomycete Phytophthora capsici is a destructive pathogen of a wide range of vegetable hosts, especially peppers and cucurbits. A 94.17-Mb genome assembly was constructed using PacBio and Illumina data and annotated with support from transcriptome sequencing (RNA-Seq) reads.

Correction: Dengue virus serotype 2 infection alters midgut and carcass gene expression in the Asian tiger mosquito, Aedes albopictus.

[This corrects the article DOI: 10.1371/journal.pone.0171345.].

Draft Genome Sequence of a Gordonia sp. Isolated from the Soil of a Red Alder Plant.

A Gordonia species was cultured from soil of a red alder (Alnus rubra) plant. Here we present the assembled and annotated genome sequence to aid investigations into the potential of this organism as a symbiont and comparative studies of the genus Gordonia.

Strategies for optimizing BioNano and Dovetail explored through a second reference quality assembly for the legume model, Medicago truncatula.

BACKGROUND: Third generation sequencing technologies, with sequencing reads in the tens- of kilo-bases, facilitate genome assembly by spanning ambiguous regions and improving continuity. This has been critical for plant genomes, which are difficult to assemble due to high repeat content, gene family expansions, segmental and tandem duplications, and polyploidy. Recently, high-throughput mapping and scaffolding strategies have further improved continuity. Together, these long-range technologies enable quality draft assemblies of complex genomes in a cost-effective and timely manner. RESULTS: Here, we present high quality genome assemblies of the model legume plant, Medicago truncatula (R108) using PacBio, Dovetail Chicago (hereafter, Dovetail) and BioNano technologies. To test these technologies for plant genome assembly, we generated five assemblies using all possible combinations and ordering of these three technologies in the R108 assembly. While the BioNano and Dovetail joins overlapped, they also showed complementary gains in continuity and join numbers. Both technologies spanned repetitive regions that PacBio alone was unable to bridge. Combining technologies, particularly Dovetail followed by BioNano, resulted in notable improvements compared to Dovetail or BioNano alone. A combination of PacBio, Dovetail, and BioNano was used to generate a high quality draft assembly of R108, a M. truncatula accession widely used in studies of functional genomics. As a test for the usefulness of the resulting genome sequence, the new R108 assembly was used to pinpoint breakpoints and characterize flanking sequence of a previously identified translocation between chromosomes 4 and 8, identifying more than 22.7 Mb of novel sequence not present in the earlier A17 reference assembly. CONCLUSIONS: Adding Dovetail followed by BioNano data yielded complementary improvements in continuity over the original PacBio assembly. This strategy proved efficient and cost-effective for developing a quality draft assembly compared to traditional reference assemblies.

Dengue virus serotype 2 infection alters midgut and carcass gene expression in the Asian tiger mosquito, Aedes albopictus.

BACKGROUND: The Asian tiger mosquito, Aedes albopictus is currently an important vector for dengue, chikungunya and Zika virus, and its role in transmission of arthropod-borne viruses (arboviruses) may increase in the future due to its ability to colonize temperate regions. In contrast to Aedes aegypti, the dominant vector of dengue, chikungunya and Zika virus, genetic responses of Ae. albopictus upon infection with an arbovirus are not well characterized. Here we present a study of the changes in transcript expression in Ae. albopictus exposed to dengue virus serotype 2 via feeding on an artificial bloodmeal. METHODOLOGY/PRINCIPAL FINDINGS: We isolated midguts and midgut-free carcasses of Ae. albopictus fed on bloodmeals containing dengue virus as well as controls fed on virus-free control meals at day 1 and day 5 post-feeding. We confirmed infection of midguts from mosquitoes sampled on day 5 post-feeding via RT-PCR. RNAseq analysis revealed dynamic modulation of the expression of several putative immunity and dengue virus-responsive genes, some of whose expression was verified by qRT-PCR. For example, a serine protease gene was up-regulated in the midgut at 1 day post infection, which may potentially enhance mosquito susceptibility to dengue infection, while 14 leucine-rich repeat genes, previously shown to be involved in mosquito antiviral defenses, were down-regulated in the carcass at 5 days post infection. The number of significantly modulated genes decreased over time in midguts and increased in carcasses. CONCLUSION/SIGNIFICANCE: Dengue virus exposure results in the modulation of genes in a time- and site-specific manner. Previous literature on the interaction between mosquitoes and mosquito-borne pathogens suggests that most of the changes that occurred in Ae. albopictus exposed to DENV would favor virus infection. Many genes identified in this study warrant further characterization to understand their role in viral manipulation of and antiviral response of Ae. albopictus.

De novo assembly of a tadpole shrimp (Triops newberryi) transcriptome and preliminary differential gene expression analysis.

Next-generation sequencing techniques, such as RNA sequencing, have provided a wealth of genomic information for nonmodel species. Transcriptomic information can be used to quantify the patterns of gene expression, which can identify how environmental differences invoke organismal stress responses and provide a gauge in predicting species adaptability. In our study, we used RNA sequencing to characterize the first transcriptome from a naupliar tadpole shrimp (Triops newberryi) to identify the genes expressed during the early life history stages and which could be important for future genomic studies. RNA was extracted from naupliar T. newberryi that were reared in a laboratory-controlled setting and in two different water types, a native and a non-native condition. A total of six replicates, three per condition, were sequenced with the Illumina Hi-Seq 2000 achieving 365 M 50-nt reads. High-quality reads were produced and de novo assembly was used to construct a T. newberryi transcriptome that was approximately 24.8 M base pairs. More than 10 000 peptides were predicted from the assembly, and genes were sorted into gene ontology categories. The use of different water conditions allowed for a preliminary differential gene expression analysis in order to compare the changes in gene expression between conditions. There were 299 differentially expressed genes between water conditions that might serve as a focal point for future genomic studies of Triops acclimation to different environments. The Triops transcriptome could serve as vital genomic information for additional studies on Branchiopod crustaceans.

A survey of the sorghum transcriptome using single-molecule long reads.

Alternative splicing and alternative polyadenylation (APA) of pre-mRNAs greatly contribute to transcriptome diversity, coding capacity of a genome and gene regulatory mechanisms in eukaryotes. Second-generation sequencing technologies have been extensively used to analyse transcriptomes. However, a major limitation of short-read data is that it is difficult to accurately predict full-length splice isoforms. Here we sequenced the sorghum transcriptome using Pacific Biosciences single-molecule real-time long-read isoform sequencing and developed a pipeline called TAPIS (Transcriptome Analysis Pipeline for Isoform Sequencing) to identify full-length splice isoforms and APA sites. Our analysis reveals transcriptome-wide full-length isoforms at an unprecedented scale with over 11,000 novel splice isoforms. Additionally, we uncover APA of ∼11,000 expressed genes and more than 2,100 novel genes. These results greatly enhance sorghum gene annotations and aid in studying gene regulation in this important bioenergy crop. The TAPIS pipeline will serve as a useful tool to analyse Iso-Seq data from any organism.

Complete Genome Sequence of Streptomyces venezuelae ATCC 15439, Producer of the Methymycin/Pikromycin Family of Macrolide Antibiotics, Using PacBio Technology.

Here, we report the complete genome sequence of Streptomyces venezuelae ATCC 15439, a producer of the methymycin/pikromycin family of macrolide antibiotics and a model host for natural product studies, obtained exclusively using PacBio sequencing technology. The 9.03-Mbp genome harbors 8,775 genes and 11 polyketide and nonribosomal peptide natural product gene clusters.

Repertoire of virus-derived small RNAs produced by mosquito and mammalian cells in response to dengue virus infection.

RNA interference (RNAi) is the major defense of many arthropods against arthropod-borne RNA viruses (arboviruses), but the role of RNAi in vertebrate immunity to arboviruses is not clear. RNA viruses can trigger RNAi in vertebrate cells, but the vertebrate interferon response may obscure this interaction. We quantified virus-derived small RNAs (vRNAs) generated by mosquito (U4.4) cells and interferon-deficient (Vero) and interferon-competent (HuH-7) mammalian cells infected with a single isolate of mosquito-borne dengue virus. Mosquito cells produced significantly more vRNAs than mammalian cells, and mosquito cell vRNAs were derived from both the positive- and negative-sense dengue genomes whereas mammalian cell vRNAs were derived primarily from positive-sense genome. Mosquito cell vRNAs were predominantly 21 nucleotides in length whereas mammalian cell vRNAs were between 12 and 36 nucleotides with a modest peak at 24 nucleotides. Hot-spots, regions of the virus genome that generated a disproportionate number of vRNAs, overlapped among the cell lines.

A transcriptomic approach to elucidate the physiological significance of human cytochrome P450 2S1 in bronchial epithelial cells.

BACKGROUND: Cytochrome P450 2S1 (CYP2S1) is an orphan P450 with an unknown biological function. Data from our laboratory and others suggest that CYP2S1 may have an important physiological role in modulating the synthesis and metabolism of bioactive lipids including prostaglandins and retinoids. CYP2S1 expression is elevated in multiple epithelial-derived cancers as well as in the chronic hyperproliferative disease psoriasis. Whether CYP2S1 expression in proliferative disease is protective, detrimental, or neutral to disease progression remains to be determined. Two human bronchial epithelial cells (BEAS-2B) were constructed to represent chronic depletion of CYP2S1 using short-hairpin RNA (shRNA) silencing directed toward the 3'UTR (759) and exon 3 (984) of the CYP2S1 gene and compared with a non-targeting shRNA control (SCRAM). Both CYP2S1 mRNA and protein were depleted by approximately 75% in stable cell lines derived from both targeted shRNA constructs (759 and 984). To elucidate the biological significance of CYP2S1, we analyzed transcriptome alterations in response to CYP2S1 depletion in human lung cells. RESULTS: RNA-sequencing (RNA-seq) analysis was performed to compare the transcriptome of the control (SCRAM) and the CYP2S1-depleted (759) BEAS-2B cell lines. Transcriptomes of the replicates from the two cell lines were found to be distinct populations as determined using Principal Component Analysis and hierarchical clustering. Approximately 1000 genes were differentially expressed in response to CYP2S1 depletion. Consistent with our previous phenotypes, DAVID analysis revealed altered regulation in key pathways implicated in cell proliferation and migration. Transcriptomic profiles were also consistent with the metabolism of proposed endogenous substrates. Pathway analysis also revealed significant expression changes within mTOR signaling, a critical pathway in cell growth. To determine whether these changes manifest as altered cell size, cell diameter and volume were calculated, revealing that CYP2S1 depletion promotes cell growth in BEAS-2B cells. CONCLUSIONS: These data suggest that pathway analysis of sequence-based gene expression is a powerful method to identify pathways and phenotypic alterations in response to changes in orphan enzyme expression. Our results suggest a novel role for CYP2S1-mediated metabolism in modulating BEAS-2B cell size. These findings warrant further studies on CYP2S1 regulated pathways to elucidate potential substrates of CYP2S1.

Transcriptome profiling of cytokinin and auxin regulation in tomato root.

Tomato is a model and economically important crop plant with little information available about gene expression in roots. Currently, there have only been a few studies that examine hormonal responses in tomato roots and none at a genome-wide level. This study examined the transcriptome atlas of tomato root regions (root tip, lateral roots, and whole roots) and the transcriptional regulation of each root region in response to the plant hormones cytokinin and auxin using Illumina RNA sequencing. More than 165 million 1×54 base pair reads were mapped onto the Solanum lycopersicum reference genome and differential expression patterns in each root region in response to each hormone were assessed. Many novel cytokinin- and auxin-induced and -repressed genes were identified as significantly differentially expressed and the expression levels of several were confirmed by qPCR. A number of these regulated genes represent tomato orthologues of cytokinin- or auxin-regulated genes identified in other species, including CKXs, type-A RRs, Aux/IAAs, and ARFs. Additionally, the data confirm some of the hormone regulation studies for recently examined genes in tomato such as SlIAAs and SlGH3s. Moreover, genes expressed abundantly in each root region were identified which provide a spatial distribution of many classes of genes, including plant defence, secondary metabolite production, and general metabolism across the root. Overall this study presents the first global expression patterns of hormone-regulated transcripts in tomato roots, which will be functionally relevant for future studies directed towards tomato root growth and development.